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Abstract 


In this study, the influence of different microwave powers (90, 270, and 450 W) and pulsation periods, i.e., 
On/Off times (30/60 and 30/120 ss!) on drying rate and seed breakage of “Nemat” and “Hashemi” paddy 
varieties” was investigated. According to the results, the Midilli et al.'s model showed the best prediction 
accuracy for drying rate of the "Nemat" variety in the pulsation period of 30/60. Furthermore, the two-term 
model was found as the best model for "Hashemi" variety in the pulsation periods of 30/60 and 30/120 and for 
"Nemat" variety in the pulsation period of 30/120.The 270 W microwave power and 30/120 pulsation period can 
be recommended as a final conclusion of this study for drying “Hashemi” paddy variety. The breakage 
percentage at this condition was 19.1%. Breakage percent of the Nemat variety was more than 40% at all of the 
conducted measurements, indicating that this variety is not suitable for microwave drying. The final concluding 
massage of this study is that a pre-test should be conducted before applying the microwave for paddy drying. 
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Introduction 


Rice (Oryza sativa L.) is a valuable source 
of energy and nurtents, which is consumed by 
almost half of the world population. It has 
been producing in tropical and sub-tropical 
countries such as India, Thailand, the 
Philippines and several other countries. The 
latest data shows that world paddy production 
in 2016 was745.5 million tons (495.2 million 
tons milled rice), FAO (2016). 

In Iran, rice is an ancient crop which is 
widely grown on areas of about 615000 ha 
with an annual production of about 3 million 
tons. Main areas of rice cultivation in Iran are 
located in the Northern provinces, Guilan and 
Mazandaran, producing 75% of the total 
cultivated rice crop in this country (Alizade et 
al. 2006). 

Agricultural crops commonly contain a 
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high level of moisture and microorganisms at 
the harvest time. After harvesting, the moisture 
content of rice paddy is in the range of 25% - 
28% (wet basis) and even higher during the 
rainy season. For this reason, immediate 
drying is a requirement in _ postharvest 
processing to avoid quality losses of these 
perishable agricultural products (Balbay et al. 
2012, Al-Harahsheh et al. 2009, Soysal 2004). 
Drying is one of the most widespread methods 
for post-harvest protection of agricultural 
products such as paddies for allowing quick 
preservation (Dadali et al. 2008, Doymaz and 
Kocaygit 2011, Discala et al. 2013). Rice 
paddy is typically dried to reduce the moisture 
content to 11% or lower for a safe storage 
before a milling process. However, if the 
moisture content in paddy is too low, the grain 
will be brittle in the milling process. This can 
lead to a higher fraction of broken kernels. 
Keeping the rice paddy at optimized moisture 
content can prolong storage time and prevent 
mould growth (Cheenkachorn 2007). 

Several drying approaches are used in the 
drying of foodstuff. The use of microwave 
technique in the drying of products has 
recently become common due to the quick and 
effective heat distribution in the product (Li et 
al. 2009, Alibas 2010, Dong et al. 2011). In 
this regard, many mathematical models have 
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been formulated in order to describe the thin 
layer drying process (Kardum et al. 2011). 
Thin layer models can be categorized as 
theoretical, semi-empirical and empirical 
models (McMinn 2006, Alibas 2014). 

The aims of this study were: I) to 
investigate the kinetics of thin layer 
microwave drying of two more conventional 
paddy varieties in Iran, ie., Nemat and 
Hashemi, II) derivation of the semi-empirical 
models for the moisture ratio and seed 
breakage percentage as a function of drying 
power and time, III) evaluation of the 
conventional thin layer drying models in 
application of the microwave drying 
conditions, which have not been examined 
before. 


Materials and Methods 


The paddy grain used in this study was 
obtained from Rice Research Institute located 
in the Sari Agricultural Sciences and Natural 
Resources University, Iran. The evaluated 
paddy varieties (Nemat and Hashemi) are 
common varieties of paddy in the north of Iran 
(Zareiforoush et al. 2009). Initial moisture 
content (MC) of the paddy was determined by 
drying the samples in a laboratory oven at 
103°C for 48 hr (Sacilik et al. 2003). Initial 


moisture content of "Nemat" and "Hashemi" 
were 20.48% and 24.5% w.b, respectively. 
Microwave drying experiments —_ were 
performed in a domestic digital microwave 
oven (MW-F304ADY-S.Media, China). The 
microwave oven had 5 different microwave 
generation power stages between 90 and 900 
W. The inside area of the device was 520mmx 
467mmx 335mm with a rotation plate of a 
300mm in diameter at the base of the oven. 
Microwave drying tests were carried out at 
three different microwave generation powers 
of 90, 270, and 450 W. The mass of each 
paddy sample was 30 g. Moisture loss was 
determined by weighting the samples using a 
digital balance with 0.01 precision. For a 
pulsed mode, the magnetron was alternately 
turned on and off corresponding to the 
specified pulsation periods. Paddy was hulled 
with a laboratory husker. Rice kernels with 


lengths smaller than 75% of that of whole 
kernels were chosen as broken rice (Kalantari 
and Eshtevad, 2013). The whole experimental 
period including several 30/60 or 30/120 ss"! 
On/Off pulsation periods was 450 s in all of 
the conducted experiments in this study. 


Mathematical formulation 


Moisture ratio of the samples after drying 
was calculated using the following equation 
(Manikantun et al. 2012): 

(M, = M,) 

where MR is dimensionless moisture ratio, 
M is moisture content at time t, Mo is initial 
moisture content, and Me is equilibrium 
moisture content on the wet basis. 

In this study, the regression coefficient (R7) 
was considered as the main criterion for 
choosing the most appropriate model for 
describing the microwave drying curves. This 
correlation factor can be used to examine the 
linear relationship between the measured and 
the estimated values, which can be calculated 
using the following expression (Alibas 2014): 
Tiee( Me, i Moment? ~ (Mp = Mg.) 


pe — 
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where R? is the coefficient of correlation, 
M ae : ; 
Rexp.i is experimental moisture ratio obtained 


by the measurements, Mae is estimated 
moisture ratio, and N is total number of 
observations. 

The standard error of estimate (SEE) 
indicates information on the long-term 
performance of the actual deviation between 
predicted and measured values term by term. 
The ideal value of SSE is "zero". The SEE is 
given as Equation 3 (Alibas 2014): 

[Pes (Marans. Mr...) (3) 


SEE = 
N- Ty 


Where n; is the number of constants. 


(2) 
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Root mean square error (RMSE) provides 
information on the short-term performance 
which can be computed from the following 
equation (Alibas 2014): 


RMSE = |e) Fi (Mapa (4) 


To choose an appropriate model for 
describing the drying kinetics of the examined 
paddies, ten empirical and semi-empirical thin- 
layer drying models were evaluated (Table 1). 
However, some of the proposed models 
presented in this table have been constructed 
based on the previous models with some minor 
modifications .In this research, MATLAB 
2013 has been used for model analyzing. 


Result and Discussion 


Effect of different powers and drying time on 
moisture ratio 

The simultaneous influence of different 
powers and drying time on the moisture ratio 
(MR) in two verities of "Nemat" and 
"Hashemi" and two pulsation periods of 30/60 
and 30/120 are shown in Fig.1 (a) to (d). 
According to the results presented in Fig.1 (a) 
to (d), the moisture ratio was decreased 
simultaneously with the power increment and 
drying time. Initially, moisture ratio was 
decreased with a high rate due to the high 
moisture content of the paddy. Meanwhile, 
moisture ratio in the 30/60 pulsation period 
was decreased with a larger slope in 
comparison with that of 30/120 pulsation 
period. A prior investigation conducted by 
Cheepsathit and Pattala (2005) indicated that 
the product temperature at the pulsation period 
of 30/60 was higher than that of 30/120. 
Higher microwave powers led to _ higher 
temperature gradients inside the grain allowing 
water to faster evaporate during the power-on 
time, also increasing diffusion of water during 
the power-off time. In contrast, longer off 
times resulted in temperature decrease due to 
evaporative cooling, convection and radiation 
heat losses from paddy external surfaces. 
These latter results obtained in this study were 
similar to the previous observations on paddy 


drying (Yongsawasdigul and Gunasekaran 
1996, Cheenkachorn 2007). 


(d) 
Fig.1. Effect of different powers and drying time on 
moisture ratio in ““Nemat” and “Hashemi”’ variety at 
30/60(ss*!) and 30/120(ss"') pulsation periods, (a); (b); 
(c); (d), respectively. 


Evaluating the existing drying models 
Evaluation of the 10 different well-known 
thin-layer drying models (summarized in 
Table 1) is presented in Table 2and 3 for 
“Nemat” and “Hashemi” varieties respectively 
at three different microwave powers. In these 
tables, the values of the standard error of 
estimate (SEE), coefficient of correlation (R’), 
and root mean square error (RMSE) are 
presented for two examined thin-layer drying 
models (i.e., Midili et al.'s model and two-term 


Table 2. Statistical results obtained from thin-layer drying 
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model). The best drying model was selected in 
this study based on the R?, RMSE and SEE 
values for the microwave drying powers of 90, 


270, and 450 W and for the 30/60 and 
30/120pulsation periods (Doungporn et al., 
2012). 


Table 1. The most conventional mathematical thin-layer drying models. 


Model No. Model name Model equation Eq No. 
1 Lewis (Doymaz and Ismail 2011) MR = exp(—kt) (5) 
2 Page (Jangam et al. 2008) ME = exp(—kt") (6) 
3 Henderson and Pabis (Pehlivan and Togrul 2004) MR = a exp(—kti (7) 
4 Logarithmic (Kingsly et al. 2007) MR = exp(—kt}+c (8) 
5 Two-term (Demirhan and Ozbek 2011) MR = aexp(—kpt!) + bexp(—k,t} (9) 
6 Wang and Singh (Demirhan and Ozbek 2011) MR = 1+ at + bt? (10) 
7 Midilli et al. (Midilli et al. 2002) MR = aexp(—kt")}+bt (11) 
8 Weibull distribution (Babalis 2006) MR = a— bexp[—(ht")] (12) 
9 Aghlasho et al. (Aghlasho et al. 2009) MR = exp( =) (13) 

ee ‘ = Ay 
10 Logistic (Alibas 2014) MR = (tacit) (14) 


MR is moisture ratio, k, ko, ki, kz, a, ao, b, c are drying constant and t is drying time. 


Regarding the results summarized in Table 
2 for (R’), (RMSE) and (SEE) values, the 
model proposed by Midilli et al had a better 
prediction for the Nemat variety at 30/60 
pulsation period, whereas the two-term model 
indicated a better estimation at 30/120 
pulsation period. The average of the standard 
error of estimate (SEE), the coefficient of 
correlation (R’), and root mean square error 
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0.01481 0.08605 0.8894 


0.9119 


(RMSE) for the model of Midilli et al at 30/60 
pulsation period were equal to 0.00053854, 
0.9840, and 0.013130667, respectively. For 
two-term model at 30/120 pulsation period, the 
average of the standard error of estimate 
(SEE), the coefficient of correlation (R?), and 
root mean square error (RMSE) were equal to 


0.000940283, 0.9783, and 0.03847, 
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information given in Table 3, the average of 
standard error of estimate (SEE), coefficient of 
correlation (R?), and root mean square error 
(RMSE) for the “ Two-term model” at 30/60 
pulsation period were equal to 0.000553233, 
0.9726, and 0.014376333, respectively. For 
two-term model at 30/120pulsation period, the 
average of the standard error of estimate 
(SEE), the coefficient of correlation (R7), and 
root mean square error (RMSE) were equal to 
0.000145233, 0.9974, and 0.014342333, 
respectively. 

Drying constant and coefficient values for 
the Midili et al.'s and “T'wo-term model” have 
been computed using the Matlab software. The 
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obtained data for the parameters existing in 
their models including k, ko, ki, n, a, and b are 
given in Table 4 for the Nemat and Hashemi 
varieties. 

Variation of the moisture ratio (MR) as a 
function of time (t) was given in Fig.2 (a) to 
(d) for two examined paddy grains. Prediction 
of the Midili et al.'s model for “Nemat” variety 
at 30/60(ss"!) pulsation periods is presented in 
Fig 2(a). The results presented in this figure 
indicate a good prediction of this model at the 
lower power-off time, whereas the ““Two-term 
model” proposed by Demirhan and Ozbek 
(2011) had a better prediction for the longer 
off-times, i.e., 120 s off-time. It is interesting 
that the “Tow-term model” showed the best 
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prediction at 120 s off-time for all of the 
microwave powers for both "Nemat" and 
"Hashemi" verities. The latter results can be 
observed in Fig.2 (b) to (d). The average 
coefficient of correlation (R*) for the “Two- 
term model” was the highest amount within 
the 10 examined models, (Table 2). For 
“Two-term model’, the value of the coefficient 
of correlation (R?) was very close to "1", 
meaning that the estimated data corresponded 
well with the experimental data. 

Considering the results presented in 
Fig.2(a) to (d), moisture ratio decreased with 
increasing the power and drying time. This 
result has a good agreement with the previous 
observations on drying of paddy, 
(Yongsawasdigul and Gunasekaran 1996, 
Cheenkachorn 2007,Cheepsathit and Pattala 
2005,Zarein et al., 2013 and Chungcharoen et 


al., 2015). 

Based on the results obtained in this study, 
moisture ratio in the 30/60 pulsation period 
was decreased with a higher rate in 
comparison with the 30/120 pulsation period 
only at higher power rates, ie., 450 W. At 
lower powers rate (90 and 270 W), influence 
of the 30/60 and 30/120 pulsation periods was 
not significant, (Table 5). This result is also 
presented in Fig 3(a) to (d), indicating that 
moisture transfer inside the paddy occurs 
faster at the higher microwave powers and 
lower microwave power off-times. This 
phenomenon can be explained by the more 
continuously heat generation within the paddy 
at higher microwave powers, yielding a larger 
vapor pressure difference between the center 
and surface of the grain. 


Table 4. Statistical results and coefficients obtained from thin-layer drying models for the different microwave power and 
pulsation period for ““Nemat” and ‘“Hashemi’’ varieties. 


Pulsation Power 
Model __ period (Watt) SEE RMSE R? a b K n Ko Ki 
(s/s) 
30/60 0.0001161 0.007619 0.9817 1.001 0.0007988 0.0007667 1.089 : : 
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3 30/120 - - - : : - : E E 
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z 30/120 , : : . : . . : : 
a 30/60 i : : ; . : : . 2 : 
: 30/120 0.001454 0.0001454 0.9556 0.0161 0.9839 : : 0.1186  0.0001408 
s 30/60 : : . . : : : ‘ . 
o ¢ 
4 5 270 : 
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A 16 
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7 
= 30/60 90  0.0001536 0.008762 0.9477 -1.484. 2.48 : -  0,0006935 0.004564 
2g 30/120 0.0001547 0.008774 0.9999 0.7071 0.2929 : - 0.000907 0.001051 
: EZ 30/60 -270-0.001373 0.02621 0.9706 ocgyy (1.057 : - 0.007436 0.007637 
§ $ — 30/120 0.0001764 0.0261 ~—«-0.9935 ‘1.038 __—-0.03759 , -  0,0006226 0.1161 
Z a 
Z 30/60 450 0.001331 0.008157 0.9994. 3.372 «4.373 ; -  1.447e-07 0.003113 
= 30/120 0.001046 0.008153 0.9988 1.055 __-0.05509 ; - 0.001236 0.112 
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Fig.2. Moisture ratio versus drying time, comparing prediction curves (~__) with the experimental data ( single 

points): a) Model of Midilli et al for the ““Nemat” variety and 30/60(ss"') pulsation periods; b) Two-term model 

for the “Nemat” variety and 30/120(ss‘!) pulsation periods; c) Two-term model for the “Hashemi” variety and 
30/60(ss"!) pulsation periods; and d) Two-term model for the “Hashemi” variety and 30/120(ss"') pulsation 


periods. 
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Table 5. Moisture ratio at the end of drying time for all powers and variety 
Power (Watt) 
Variety Pulsation period (ss') 90 270 450 


Weacai 30/60 0.91 0.72 0.36 
30/120 0.92. 0.72 0.49 

Hashemi 30/60 0.93 0.76 0.73 
30/120 0.94 0.79 0.61 


Effect of different powers and drying time on 
seed breakage percentage 

Simultaneous influence of different powers 
and drying time on the seed breakage in two 
"Nemat" and "Hashemi" varieties and two 
pulsation periods of 30/60 and 30/120 are 
shown in Fig.3 (a) to (d). 

According to Fig.3 (a) to (d), the seed 
breakage was increased simultaneously with 
increasing the microwave power and drying 
time. The experimental data showed that both 
microwave power and pulsation period affect 
the broken fractions. For the pulsation period 
of 30/60, the higher broken fraction was 
observed at higher microwave powers, e.g., 
450 W. These results might due to the higher 
energy absorption of the paddy grains which 
led to a rapid increase in the grain temperature. 
However, a rapid increase of grain temperature 
resulted in the higher rate of water removal. 
Excessively high moisture content removal 
most likely causes a high fraction of broken 
kernels. The results were similar to those of 
previous studies conducted by Soponronarit et 
al. (1996), Taweerattanpanich et al. (1999) and () 
Cheenkachorn (2007). The obtained empirical —— 
expressions for the seed breakage percentage 
(SB %) as a function of microwave power and 
drying time together with the moisture ratio 
(MR) are given in Table 6. These expressions 
have been derived using MATLAB software 


Conclusion ; : (d) . 
In this investigation, the influence of Fig.3.Effect of different powers gud dayne time on 
. . ‘ seed breakage percent in ‘“Nemat” variety and 
different microwave powers and pulsation “Hashemi” variety at 30/60(ss) and 30/120(ss"!) 
periods on the moisture ratio and _ seed pulsation periods, (a); (b); (c); (d), respectively. 


breakage of two paddy _ varieties was 
investigated. According to the obtained results 
in this study, the following conclusions can be 
summarized. 
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The minimum moisture ratio and _ the 
maximum seed breakage percentage were 
observed at the microwave power of 450 
W and pulsation period of 30/60. This can 
be explained by the higher heat generation 
inside the grain and consequently higher 
evaporation rate at the grain surfaces, 
yielding larger pressures and thermal 
gradients inside the grain. 

The maximum moisture ratio and _ the 
minimum seed breakage percentage were 
obtained at the microwave power of 90 W 
and pulsation period of 30/120. 

At all of the examined microwave powers 
(90, 270 and 450 W), the “Midilli et al.'s 
model” and “T'wo-term model” have been 
found to be the best model for the 
"Nemat" variety in the 30/60 and 30/120 
pulsation periods, respectively. 

For all of the examined microwave 
powers, the “Two-term model” has been 
found to be the best model for the 
"Hashemi" variety in the 30/60 and 
30/120 pulsation periods. This model has 
been selected using the _ statistical 
information including R*, SEE, and 
RMSE and_ simultaneously comparing 
with the experimental data. 

Moisture ratio in the 30/60 pulsation 
period was decreased with a higher rate in 
comparison with 30/120 pulsation period 
at higher microwave power, i.e., 450 W 
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for both examined varieties. This result is 
in agreement with the prior investigation 
carried out by Cheepsathit and Pattala 
(2005). 

The Nemat variety is not recommended 
for a microwave drying based on the 
results obtained in this study for all of the 
examined microwave powers and _ test 
conditions. The minimum breakage 
percentage was obtained about 40% at 90 
W and 270 W microwave powers and 
30/120 pulsation period. 

For the Hashemi variety, the minimum 
breakage percentage was obtained equal 
to 17.6% at 90 W and 30/120 pulsation 
period. The breakage percentage at 270 W 
and 30/60 and 30/120 pulsation periods 
were obtained equal to 20% and 19.1%, 
respectively. Therefore overall and 
simultaneously considering the breakage 
percentage and the drying rate (Figs.3 (a) 
to (d),the 270 W microwave power and 
30/120 pulsation period can _ be 
recommended as a final concluding results 
of this study. The 270 W microwave 
power with 30/60 pulsation period is not 
recommended due to the larger energy 
consumption and no significant drying 
rate increment and breakage percenaget 
decrement in comparison with 270 W and 
30/120 pulsation period. 
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